arXiv:1505.06121v2 [hep-ph] 23 Jul 2015 


IFIC/15-30 


(S'[/(5)-inspired double beta decay 


Renato M. FonsecaE 


and Martin 


Hirscliil 


AHEP Group, Instituto de Fisica Corpuscular - C.S.I.C./Universitat de Valencia 
Edificio de Institutos de Paterna, Apartado 22085, E-46071 Valencia, Spain 

Abstract 


The short-range part of the neutrinoless double beta amplitude is generated via the exchange of 
exotic particles, such as charged scalars, leptoquarks and/or diquarks. In order to give a sizeable 
contribution to the total decay rate, the masses of these exotics should be of the order of (at most) 
a few TeV. Here, we argue that these exotics could be the “light” (i.e weak-scale) remnants of 
some B — L violating variants of SU{5). We show that unification of the standard model gauge 
couplings, consistent with proton decay limits, can be achieved in such a setup without the need 
to introduce supersymmetry. Since these non-minimal S'C/(5)-inspired models violate B — L, they 
generate Majorana neutrino masses and therefore make it possible to explain neutrino oscillation 
data. The “light” coloured particles of these models can potentially be observed at the LHC, and 
it might be possible to probe the origin of the neutrino masses with AL = 2 violating signals. 
As particular realizations of this idea, we present two models, one for each of the two possible 
tree-level topologies of neutrinoless double beta decay. 
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I. INTRODUCTION 


Neutrinoless double beta decay (0i//3/9) is usually considered as a probe to constrain (or 
measure) the Majorana neutrino mass scale. However, a non-zero amplitude for decay 
is generated in any extension of the standard model with lepton number violation (for a 


recent review see, for example [ 1 |) and, in general, the mass mechanism is not necessarily 
the dominant contribution to the decay rate. 

The short-range part of the 0uPl3 decay rate is generated via the exchange of exotic 
particles, such as leptoquarks and/or diquarks plus possibly some exotic fermions. The 
“heavy” mediators in these diagrams must have masses of at most a few TeV in order 
to give a sizeable contribution to the decay rate, and thus they can potentially be 


produced and studied at the LHC [3|, . A complete list of all (scalar) mediated short-range 

contributions to decay has been given in j^. One can understand the results of this 

work as a bottom-up reconstruction of all possible particle physics models with tree-level 
00(3(3 decay. In this paper we take the opposite approach and study instead the 00(3(3 decay 
in a top-down approach; in other words, we explore the possibility that some of these exotics 
could be the “light” remnants of multiplets in some unihed theory. 

Minimal SU{5) accidentally conserves H — L j^. Once SU{5) is broken to the standard 
model group (Gsm), baryon number violating processes such as proton decay occur, but the 
combination B — L is still conserved. Thus, neutrinos are as massless in minimal SU{5) 
as in the standard model (SM) and remain so also after SU{5) breaking. Adding SU{5) 
singlets a type-I seesaw mechanism j7-10| can be generated. While certainly theoretically 
attractive to many, this scenario leaves as its only prediction that neutrinoless double beta 
decay should be observed. ^ 

Allowing for larger SU{5) multiplets, however, B — L is no longer conserved in SU{5). 
Perhaps adding a scalar 15 is the simplest way to obtain a model with B — L violation in 
the SU{5) symmetric phase . Indeed, it is possible to write down the Yukawa interaction 
Si. ■ 15 ■ 5^ as well as a scalar trilinear term 5 ■ 15* ■ 5, with the 5i. containing the SM 

and L fermion helds and the 5 containing the SM Higgs doublet H . In terms of Gsm 
representations, we have both the LLS' 1^34 interactions, where C 15 

therefore a type-H seesaw lll-ll4| will be generated, yielding the effective operator LLHH 
(or Si. • Si? • 5 ■ 5 in terms of SU{5) helds). Many other examples oi B — L violating SU{5) 
models can be constructed using larger representations. 

Extrapolating the SM gauge couplings to larger energies with only the SM particle content 
fails to achieve unihcation. Thus, no consistent model of grand unihcation can be built 
without adding new particles and/or interactions at some so-far unexplored energy scale. 


^ Albeit without fixing the 0vj3(3 decay half-life, unless further assumptions about the active neutrino mass 
spectrum are made. 

^ We use S to denote a scalar and ip for a fermion, subscripts are the transformation properties (or charge) 
under the SM group in the order SU{3)c x SU(2)l x U{1)y- For SU{5) multiplets we add a subscript 
“T” for fermions, no subscript for scalars. 
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One of the most cited possibility to achieve gauge coupling unihcation (GCU) is to extend 


the SM to the minimal supersymmetric standard model 15|]. However, it has been known for 


a long time that also many non-SUSY scenarios can lead to GCU (for an early reference with 

for an early reference with a left-right symmetric intermediate 


16 


new states at TeV, see 

stage at (10^° — GeV see [l7|). For a discussion of GCU in non-SUSY S'f/(5)-based 




models see also 18, 19 


In this paper we show that non-minimal S'U(5)-based extensions of the SM can lead to 
good GCU, provided some beyond-SM (coloured) multiplets are accidentally light. The very 
same TeV-scale remnants of these non-minimal models, responsible for GCU, could be (some 
of) the mediators of the short-range double beta decay amplitude, leading to lepton number 
violation (LNV) at the electro-weak scale. We call this “S'U(5)-inspired” double beta decay. 
Different models realizing this idea can be constructed. Since short-range diagrams for 0z//9/3 
decay always fall into one of only two possible tree-level topologies j^, we will discuss two 
particular variants, namely one model for topology-I (T-I) and one for topology-II (T-II). 

The main motivation to study these S'U(5)-inspired models is that they are experimentally 
falsifiable at the LHC and, possibly, in upcoming lepton flavour violation searches in the 
following sense. First, current limits on the half-lives of 0z//3/3 decay in ™Ge 20| and ^^®Xe 


21 


23| are of the order of (1 — 2) x 10^^ years, resulting in lower limits on the effective 


mass scale of the underlying operator in the range of ~ (2 — 2.b)g‘l^ TeV for our 
two example models, while an observation with a half-life below roughly 10^^ years implies 
Mgff < (3 — 3.8)g^g TeV. Here, ges is the mean of the couplings entering in the decay 

diagram(s), see below, with g^g required to be roughly of the order O{0.1 — 1) to give an 
observable decay rate. Negative searches at the LHC in run-H will allow the covering of 
this range of masses [i, . Second, oscillation experiments have shown^ that lepton flavour 

is violated (LFV). If our models are to explain neutrino data consistently, LFV violating 
entries in the Yukawa couplings are therefore required. These will lead to non-zero rates 
in processes such as u —?• ey, with expected rates that could be of the order of the current 
experimental limit 


25 


We should also add a disclaimer about naturalness here. Standard SU{5) suffers from 
what is known as the doublet-triplet splitting problem, i.e. the fact that the SM Higgs has 
a mass of ~ 125 GeV 26|, |27|] while the coloured triplet in the 5 must have a GUT scale 
mass. Although several solutions to this problem have been suggested (for a short review 


see [28|]), we do not concern ourselves with any particular one. Instead, we view this simply 
as a fine-tuning problem and, in fact, to make the exotic particles in our models light will 
require, in general, additional hne-tunings. We assured ourselves, however, that this can be 
done consistently for all the light states (details can be found in appendix [B]) . 

The rest of this paper is organized as follows. In section [TTl we discuss our model based 
on topology-I. It has three new multiplets when compared to minimal SU{5), one of which 
could contain a good candidate for the cold dark matter in the Universe. The model also 


See for example the recent fit 24| . 
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has a coloured octet fermion and a scalar leptoquark at the TeV scale, both of which can 
be produced at the LHC. Neutrino masses are dominated by 2-loop diagrams. Section HTTI 
discusses the model based on T-II. More light states exist in this variant model, leading to a 
more diverse phenomenology. As in model T-I, in model T-II neutrino masses are generated 
at 2-Ioop order. There are new coloured sextets, which have particularly large LHC cross 
sections and there is also a doubly charged scalar. We then close the paper with a short 
conclusion. Many of the technical aspects of our work are deferred to the appendices lAl 
IbI and [cl where we give details of the Lagrangians of the two models, discuss briefly some 
aspects of hne-tuning in SU{5) and present a table with the decompositions of larger SU{5) 
multiplets. 


II. A SIMPLE MODEL WITH A SCALAR QUADRUPLET 


We start this section with some preliminary comments. Minimal SU{5) puts the three 
standard model families of fermions into three copies of Ep and lOi;’. To generate the SM 
Dirac masses at least one scalar 5 is needed. Although not strictly speaking “minimal” we 
allow, in principle, also for the presence of a scalar 45 at the GUT scale. This 45 is added 
for the sole purpose of generating a Georgi-Jarlskog factor for fermion masses 


29 


- see 

appendix [Bl In addition, at least one scalar 24 is needed to break SU (5) to the SM group. 
As mentioned already in the introduction, the Lagrangian of this minimal model conserves 

b-lQ. 

Both our models need to introduce larger representations. For completeness we give in 
the appendix O the decomposition into SM group representations of all SU (5) multiplets 
up to the 75. Once these larger multiplets are added, B — L will be violated, as discussed 
with the example of the 15 in the introduction. 

Let us now turn to our hrst example model, which we call T-I in the following. To 
generate a topology I 0z//5/5 decay, an exotic fermion is needed. We thus add to minimal 
SU{5) three multiplets: two scalars (15 and 70) and one fermion (24p). The 15 and 24^ 
will generate Qul3l3 decay and neutrino masses, while the 70 can play the role of a dark 
matter candidate, as explained later on. Details of the Lagrangian of the model are given 
in appendix [Bl Here we only discuss the most relevant terms. Gonsider hrst the interaction 
among the 15, the 24f and the SM fermions: 




, 24 ^ 15 *+ ••• 


If we add to eq. CO a Majorana mass term for the 24^, 

= ?t7,2424j?24^ , 

an effective d = 9 operator is generated, which in SU{5) language reads 

Og (X (5ir5i?)(10F)(10i7’)(5i?5F) 


( 1 ) 


( 2 ) 


(3) 
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FIG. 1: Double beta decay short-range diagram in 5t/(5) language (left) and in the Gsm language 
(right). The fields 5 ^ and lOj? contain the standard model fermions, while T = Ss 2,1/6 and 
O = 'i/’8,i,o are the “light” pieces coming from the 15 and 2Af-, respectively. 


— see the left-hand side of fig. o. Under the SM group, the 15 and the 24^7’ break as 
5 ' 3 , 2 , i /6 + ■ ■ ■ = T + ■ ■ ■ and "08, 1,0 + ■ ■ ■ = O + ■ ■ ■, respectively. The above Lagrangian 
equations ([T]) and (E]) then contain the following terms: 

C = ylfl.d’jT + yf^Q^OT* + moOO , (4) 

with y^^"' = y^^\ y^^^ = and mo = ^24 in the SU{5) symmetric phase. 

Eq. dH) will produce an operator that generates a contribution to double beta decay via 

Otr ^ {Ld^){Q){Q){Ld^). (5) 


Here, the subscript “11” indicates the number of this AL = 2 operator in the list defined 
by Babu & Leung [30[, while the superscript “4-i” identihes the double beta decay decom¬ 
position according to the list of j^. The diagram which generates this operator is shown in 
fig- (HD, on the right. In order for this diagram to give a sizeable contribution to the total 
Offd/d decay amplitude, T and O must have masses of the TeV order. The amplitude for 
this diagram has been calculated in j^: the limit from ^^®Xe 21-23| results in Mes > 

TeV, where Meg = 5'eff is the (geometric) mean of the four couplings entering 

the diagram. In addition to this diagram, Qvjdfd decay will receive also a contribution from 
the mass mechanism. We will compare the short-range contribution and mass mechanism 
below, when discussing neutrino masses in this model. 

We now turn to the discussion of the 70 scalar representation. The standard model 
does not contain a particle candidate for the cold dark matter (DM). It is, however, quite 
straightforward to identify the basic requirements for SM multiplets to contain viable dark 
matter candidates. As shown in table H] in the appendices, the 70 is one of the smallest 
SU{5) multiplet containing an SU{2)f quadruplet, 70 = 5'i,4,i/2 + ■ ■ ■ = K + ■ ■ ■. S'!, 4,1/2 
contains one electrically neutral state, which after the breaking of SU{2)l can play the 
role of the cold dark matter. 

Of course, only particles stable over cosmologically long times can be dark matter. In the 
SU (5) phase, model T-I allows, in principle, a quartic scalar term 5 ■ 5 ■ 5* ■ 70* which could 
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FIG. 2: Gauge coupling unification in model T-I. In this calculation the new states are assumed 
to have masses around M ~ 0{1) TeV. The running includes 2-loop /3 coefficients, but does not 
consider any GUT scale thresholds. 


induce the decay —?■ hhh. To eliminate this coupling, and all other couplings linear or 
cubic in 70, one has to postulate a Z 2 symmetry. We simply assume K to be odd under 
this Z 2 , ^ while all other particles of the model are even. We will not discuss further details 


of the phenomeno logy of as a DM candidate; these have been worked out in M|. Note 


that according to M|, the neutral member of IC can play the role of DM if m^o ~ 2.4 TeV. 


Adding complete SU{5) multiplets to the SM does not change GCU, since the (3 coef¬ 
ficients of all three gauge couplings change by the same amount. However, after SU{5) 
breaking, the masses of the different Gsm multiplets contained within each SU{5) represen¬ 
tation may be different, yielding a large number of possibilities to achieve GCU. 

The f3 coefficients for the running of the gauge couplings, ® including the contributions 
from T, O and K are 




/ 326 19 28 \ 
75 3 

-1 

5 bij 

4 94 20 

-14 

7 15 88 

\ —/ 


V 6 T T/ 


at 1-loop and 2-loop order, respectively. Fig. (|2]) shows the resulting running of the inverse 
gauge couplings as a function of energy, assuming that the masses of T, O and K are of 
the TeV order. Here, we neglect the small 2-loop contributions from Yukawa couplings for 
simplicity, and we also do not consider possible corrections from GUT-scale thresholds. 

As £g. (EJ shows, the model nicely unihes with an estimated GUT scale of me — 
10^^ GeV, which gives an estimated half-life for gauge mediated proton decay of Ti/ 2 (p) ~ 


^ A non-zero vacuum expectation value (VEV) of the would break this Z 2 spontaneously and, thus, has 
to be avoided. For m\ > 0 one expects a zero VEV to be the preferred solution of the tadpole equations. 
^ In general, ignoring the effect of Yukawa interactions, the running gauge couplings gi change with the 
logarithm of the energy scale t = \ogE as follows: dgi/dt = (47r)^ -|- bijg^g^j (47r)^. 
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Iqss years, albeit with a large uncertainty. Recall that the current constraint from Super- 
Kamiokande is ^ 10^"^ years [32 . 


However, while the model is automatically safe from the gauge mediated proton decay 
diagrams, we have to impose one more constraint on the model, due to scalar mediated 
proton decay. The triple-scalar term 5 ■ 5 ■ 15*, see hg. ([3]), leads to a vertex Ss,i-i/ 3 HT. 
After electro-weak symmetry breaking this can be considered as an effective mixing of T 
with the coloured triplet Note, however, that in the full model there is more than 

one contribution to this vertex, since the field content of the model also allows the writing 
of the interaction 5 • 5 • 15* ■ 24. Furthermore, there is a second «S' 3 ,i,-i /3 inside the 45 which 
needs to be taken into consideration. As such, if we call pi and fi 2 the effective interactions 
of these two heavy colored scalars (assumed to have a degenerate mass of the order of the 
GUT scale me) with H and T, then 

/ii = 2 cos a ^ 5 . 5 . 15 * -|- 2 cos a A 5 . 5 . 15*.24 (24) — sin a Af 5 . 45 . 15 *. 24 ]^ (24) (7) 


/i 2 = cos a 


Af 5 . 45 . 15 *. 24 ]^ (24) — 2 sin a I h. 

a=l,2 \ 


a=l,2 


45-45 15 '' 




[ 45-45 15 *- 24 ], 


(24) 


( 8 ) 


a=l,2 


where a is the angle controlling the admixture of the 5 'i^ 2 ,i /2 representations in the 5 and 
45 forming the light state H, h (A) represent the trilinear (quartic) scalar couplings of 
the SU{5) helds indicated in subscript, and the index a keeps track of the different gauge 
invariant contractions of the representations. The VEV of the 24 must be in the SM-singlet 
direction. 

Since the S's^i _i /3 helds has diquark couplings, a proton decay diagram is induced by 
this 5 * 3 , 1 ,- 1/3 -H- T mixing. Consistency with proton decay limits can be converted into an 
upper limit on the sum of the two /ij couplings, p = pi -|- /i 2 , which is, very roughly, of the 
order of 

0.l\ / 2 TeVy 

\ rriT ) ' 

Here we assume that the typical value of is 0(0.1). Much smaller values, say y[^ ~ ye, 
would avoid the hne-tuning of eq. ([9]), but at the same time they would render the 0z//9/5 
diagram in hg. ([ 1 ]) unobservable. 

We note that a very similar discussion about scalar-mediated proton decay can be found 


F 

me 


< 2.5 X 10 


-6 


(9) 


m 


18 


33| . Reference [l 8 | considers also SU{5), while |33| discusses an S'O(IO) based model. 


In 33|, however, the authors argue against the existence of light colour-triplet scalars on 
the basis of this constraint, while we accept eq. ([9]) as just one (more) hne-tuning. 

Let us now discuss neutrino masses. The 15 also contains the 5 'i^ 3 ^i, which is often 


denoted as A in neutrino physics 34,l35|. Thus, after integrating out the heavy A, ehectively 


a Weinberg operator [3^ is generated. After the breaking of SU (2)^, one thus hnds a type-H 
seesaw contribution to the neutrino mass. One can estimate this contribution to be of the 
order of 

0.3 meV. (10) 


/ \ 2 /^A 

~ yljVsM—r 


m. 
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FIG. 3: Scalar-mediated proton decay in model T-I. Mixing induced via the SM Higgs VEV between 
Si, 3 ,- 1/3 and T can generate an effective diquark coupling for T. 



FIG. 4: 2-loop neutrino mass diagram in model T-I. To the left, the SU{5) origin of the diagram; 
to the right the contributions to rriu from the light states of the model. 


For = 1 and /ia = rri/^ ~ me, with a GUT scale of, see fig. ([2]), me ^ 10^^ GeV. This 
is too small to explain solar neutrino oscillations by a factor of roughly 30. 

Note that the effective /ia contains contributions from the same couplings (5 ■ 5 • 15*, 
etc.) as discussed just above for the scalar induced proton decay. The different SU{5) terms 
contribute, however, with different Glebsch-Gordon coefficients to /ia than those appearing 
in /i, see again appendix [B] for details. Thus, /i and /ia can easily take very different values. 
Of course, if /ia were to obey a limit similar to eq. ([9]), contributions to neutrino masses 
would actually be much smaller than the numerical value in eq. ffTOj) . 

The failure to have a sufficiently large neutrino mass to explain solar/atmospheric data 
via the type-II seesaw contribution ® does not imply that the model cannot explain neutrino 
oscillation data. This is due to the two-loop contribution to the neutrino mass matrix, 
involving the leptoquark and the coloured fermion, shown in fig. (j4]) (on the right). The 
left diagram shows the SU{5) origin of this diagram, the diagram to the right is in the SM 
phase. 

Using the general analysis of 2-loop diagrams of 37|] one can calculate the contribution 


® Unless mA is tuned to be much smaller than the GUT scale, which would harm the successful GCU. 







of this diagram to the neutrino mass matrix. However, since a 2-loop model based on this 
diagram has been discussed recently in 38|, we will not repeat all the details here and show 
only a few rough estimates. For mj’ ~ mo 3> hg. (j4]) gives approximately 


Nr. 




(16vr2)2 Alnv 




( 11 ) 


Nc is a colour factor {Nc = 3 for this diagram) and logarithmic terms from the loop integrals 
have been neglected for simplicity. Alnv is the mean of the masses of tut — mo- Note that 
eq. (ITT]) produces three non-zero neutrino masses, which in leading order will be proportional 
to 2,3 oc rubrurf, m^ms, While it is possible to ht quasi-degenerate neutrino masses 
with eq. fITTD . the hierarchy in down quark masses leads us to expect that neutrino masses 
follow a normal hierarchical neutrino spectrum in model T-I. For Alnv — 1 TeV, | 2 / 3 ^^| ~ 
~ 0.06 and ~ y^^ ^ 0.26, where \y^^'^\ = one obtains ~ 

(9 X 10“^, 5 X 10“^) eV, correctly reproducing the atmospheric and solar neutrino mass 
scales (for normal hierarchy). 

With two Yukawa vectors {y^^^ and y^ 2 ) contributing (dominantly) to the flavour 
structure of {m^)ij, there are four independent ratios, which we arbitrarily choose to be 
and (|/!?/|/ 3 ?, l/SVi/a?)- With these we can easily ht the observed neu¬ 
trino angles. Many solutions exist and can be found in a simple numerical scan. We 
only quote one particular example as a proof-of-principle. Choosing the absolute val¬ 
ues of the Yukawa couplings as mentioned above, {y^/y^/y^) ~ (0.03,—1.7) and 


{y 


(4) /.,(4) ^,(4) , U ) 


12 /2/a2_i 


^22/2/: 


32 


(1,1) gives all three measured neutrino angles near their best ht 


points ( 2 ^. Nevertheless, we stress again that these numerical examples give only a rough 


estimate. 

One should also compare the size of the relative contributions to the 0i'(3(3 decay rate of 


the short-range diagram (see hg. ([T])) with the one of the neutrino mass mechanism 39 


The current limits on OuPP decay from ^®Ge and ^^^Xe correspond to an upper limit on the 
ehective neutrino mass, {m„), of roughly (m^) < (0.2 — 0.4) eV, depending on the choice of 


nuclear matrix elements 40 


42|. Experiments on these nuclei with a sensitivity of 10^' years 


would probe {mu) ~ 50 meV. For a normal hierarchical neutrino mass spectrum, as in our 
example ht discussed above, much larger half-lives are expected. The short range diagram 
gives a half-life of roughly 10^^ years for TeV, corresponding to Meg >1.7 TeV 

for = \ly^u.yT' — 0-35. The latter can be completely covered by LHC searches for the 
signal /■•'/+ + as discussed below. 

We have also estimated Br(/i —>■ ey) for this model. Using the formulas 
can roughly estimate Br(/i —)■ e 7 )~ 3 x 10“^^||/|3^|/23^*p for a leptoquark mass of Mq — 1 
TeV. This is roughly of the order of the limit expected for the second phase of the MEG 


38, 39 


we 


experiment 
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A much more detailed discussion about LFV can be found in 


38 


We now discuss briehy LHC searches for model T-I. Leptoquarks have been searched at 


the LHC in run-I by both the ATLAS 43|, l4^ and the CMS 451447l| collaborations. No 


positive observation of any leptoquark state has been reported, apart from a possible 2.6(T 
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excess near ttilq — 650 GeV found in Lower limits on LQ masses depend on the 

lepton generation they couple to, and on the branching ratios of the LQs. Lower limits are 
roughly 650 — 1000 GeV |43l447l| for searches based on some pair-produced LQs, depending 
on assumptions. We also mention the recent GMS search for singly produced LQ states. 
This search excludes single production of hrst generation LQs with masses below 1730 GeV 
(895 GeV) for leptoquark couplings equal to y = 1 {y = 0.4) j48|. Run-2 will improve these 
numbers in the very near future. 

The coloured fermionic octet, O, can be pair-produced at the LHG in gluon-gluon fusion. 


We have implemented the model in SARAH ^ 


50 


and used the Toolbox environment 


51 


to generate SPheno 52|, [SSj] and MadGraph hies. We then used MadGraph5 5^ for a quick 
estimation of the cross section for pair producing O at ^/s = 8 and y/s = 13 TeV. The 
kinematics of this decay depends on the mass hierarchy between O and our leptoquark T, 
but the decays of O will always lead to a hnal state containing two hard jets plus a charged 
lepton or neutrino (i.e. missing Et)- In total, the signal will thus consist of four hard jets 
with 2,1 or 0 charged leptons. Due to the Majorana nature of O, the ratio R of the number 
of same-sign to opposite sign charged lepton events is expected to be R = 1. 

Also, we expect that the signal will contain lepton havour violating events, i.e hnal states 
plus jets (and also events with tans). This is caused by the havour structure of the 
leptoquark couplings y[l\ which is needed to explain neutrino angles and which enter in the 
decay rate of the state O. 

While there is no dedicated search for this signal at neither GMS nor ATLAS, we can 
use the results of the GMS search for right-handed neutrinos in left-right symmetric models 
551] to estimate roughly the current sensitivity of LHG to O. The GMS collaboration uses 
the hnal state 21 -|- 2j to extract lower limits on the masses of Wr and vr, but it mentions 
explicitly that in the hnal state they allow for any number of jets larger or equal to two. 
^ Below miijj ~ 2 TeV the limits derived by GMS are background dominated. Since one 
can expect that backgrounds are smaller for signals with larger number of jets, we can thus 
use the derived upper limits on a{pp —)■ lljj) to conservatively estimate lower limits on mo- 
Unfortunately, GMS decided to show limits only for mujj > 1 TeV, with upper limits on 
a{pp —)■ eejj) and a{pp —)■ ppjj) both of the order of ~ 9 fb for this lowest mass (summed 
over both lepton charges). Gomparing with their hg. (2) 55(|, one can estimate that in the 
bins [0.8,1] ([0.6,0.8]) TeV upper limits should be roughly 14 (23) fb. Gomparing with the 
calculated a{pp —)■ OO) of 24 (64) fb for mo = 900 (800) GeV and assuming a branching 
ratio Br(00 -^11 + 4j) ~ 25 %, we estimate that the current lower limit on mo should 
be roughly in the 800 — 900 GeV ballpark. More exact numbers would require a dedicated 
search by the experimental collaborations. However, for y/s = 13 TeV, cross sections will be 
much larger and we estimate a dedicated search could be sensitive up to mo — (2.1 — 2.2) 
TeV with C ~ 300 fb“^, corresponding to (16-10) events before cuts. 


^ There is a 2.8a excess in the data at rueejj — 2.2 TeV, but CMS concludes that it is not in agreement 
with expectations for a left-right symmetric model. Our model T-I cannot explain this excess either. 
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FIG. 5: Dimension-9 operators in model T-II built with SU{5) representations (left), and with 
the corresponding light Gsm pieces (right). Here, in addition to T = 2,1/6) we are using the 

shorthand notation R = X = 3 _i/ 3 , Y = Sg 3 4/3 and Z = •S’g 4 _ 2 / 3 . These operators 

generate OufijS decay (with the blue/top fields on the right diagrams) as well as n — n oscillations 
(with the red/bottom fields) — see text for a discussion. 

III. MODEL T-II: ADDING COLOURED SEXTETS 

Topology-II for decay consists of diagrams in which three scalars are interchanged 

between the six SM fermions that make up the 0z//9/9 decay operator [^. For model T-II 
we therefore add two copies of the scalar 15 plus one scalar 50 to minimal SU{5). The 
Lagrangian of the model is given in the appendix. Here we discuss only the most relevant 
terms. In the SU (5) phase the Lagrangian includes 

jCffs = ^pF,^5F,Jl5k + 10i.,a0i.,,50 + • • • , ( 12 ) 

Csss = hSf 15*15,50 + hfl5i50*50* + • • • , (13) 

where the dots stand for additional terms. Terms in eq. ffT^ . together with the corre¬ 
sponding mass terms for 15 and 50, produce the d = 9 operators shown in hg. ([5]) on the 
left. 
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We introduce the shorthand notation 


T = 


(6,3,-i).r.(6.1.i)e50^ 


(14) 

(15) 


Differently from model T-I, we will assume that there are a total of six light states: 2 copies 
of T plus one copy of -R, X, and Z. The Lagrangian in the SM phase then contains the 
following terms: 


^FFS 

^SSS 


V^lFd‘T, + yf + y\fete‘R + y'y'QiQjX + ylfuX,Y + h.c., (16) 

h\fTiTjX + h^^'>YZZ + h^^^R*Y*Z + h^^^X*X*Z + h.c.. (17) 


With these terms the diagrams in the left of hg. ([5]) produce two diagrams each in the 
SM phase of the model. These are shown in £g. ([S]) on the right. We stress again, that 
the central couplings in these diagrams, proportional to the couplings receive 

contributions from several terms of the Lagrangian in the SU (5) phase of the model and, as 
discussed in the appendix, can therefore take quite different values, despite their common 
origin. 

Consider first the diagrams on the right-hand side of fig. ([5]) containing two external 
leptons. We can easily identify them with the T-II operators classified in as 


Oi,-^^-^ = {Ld){Ld){QQ), 


O 


T-II-3 


= (eV) 


(18) 

(19) 


Here, again the subscript “11” identihes the Babu & Leung operator 3^, whereas is the 
missing AT = 2 dimension 9 operator in this list. The superscripts indicate the 0i^/3/3 decay 
label in the classihcation scheme of j^. The latter allows us to estimate the 0z//9/3 decay rates 
induced by these two diagrams. By coincidence the two operators give very similar limits: 
using again the experimental limits from |211-123| results in Meg > 2.2gJ^ TeV for each of 
the two operators. As before, Meg stands for the mean of the masses entering the diagram 
and for the definition of the mean coupling we have made and dimensionless by 
introducing g^^'^ = We stress that, in contrast to model T-I, here all light states 

appear in at least one tree-level double beta decay diagram. However, this variant model 
does not include a dark matter candidate. 

The two diagrams on the right-hand side of hg. ([5]) also produce the operators: 


Oab=2,al=o = (/ d ')(//)( ll " ll "), 
Oab=2Al=o = {QQ){d"I){QQ ). 


( 20 ) 
( 21 ) 

Both of these operators induce n — n oscillations — see for example the review [^. The 
current limit on the life-time of neutron-antineutron oscillations is Tn^n > 0.86 x 10® s M 


and, as discussed in that paper, one expects an improvement on this number of about two 
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FIG. 6: Gauge coupling unification in model T-II. In this example calculation the new states are 
assumed to have masses around M ~ 0{1) TeV. As in fig. ([2]), the running includes 2-loop (3 
coefficients, but does not consider any GUT scale thresholds. 


orders of magnitude in the future. We can estimate the life-time of neutron-antineutrino 
oscillations, Tn-n adapting the formulas from 33(] for our model. Introducing the notation 
9 = one gets from the current experimental limit a constraint 


< 0.6 X 10"® f ^ 
9 


M 


eff 


3 TeV 


GeV 


( 22 ) 


for the bottom right diagram in £g. ([5]) (in this case, M^s = (m^m|)^/®). The same 

constraint applies to from the top right diagram of hg. ([5]). This hne-tuning could only 

which 


be avoided if Z takes a mass of the order of the GUT scale, as in the model of 33 
would, however, destroy gauge coupling unihcation in model T-II discussed next. 
For model T-II we hnd the following (3 coefficients: 


h = 


/ ^ \ 


/1691 

129 

1988 \ 

30 


30 

10 

15 

11 

6 

1 ^ij ~ 

43 

10 

785 

6 

188 

13 


497 

141 

541 



\ 30 

2 

— J 


(23) 


Fig. ([6]) shows the corresponding running of the inverse gauge couplings for masses of the 
new particles M ~ 0{1) TeV. Despite the very different particle content of model T-II, 
compared to model T-I, GGU works very nicely. The estimated GUT scale in this model 
is found to be roughly ~ 5 x 10^^ GeV, lower than in model T-I. Also ac is predicted 
to be much larger than in model T-I. This can be easily understood from the f3 coefficients, 
which are larger due to the presence of colour sextets in the running. 

The best-£t point for the proton decay half-life is found to be rp^^ 0 e+ ~ 5 x 10^^ years. 


slightly below the current experimental limit 


32 


However, it is well-known that predictions 


of proton decay half-lives have a large error bar. We have repeated the exercise of estimating 
the uncertainty in the GUT scale due to the neglected GUT-scale thresholds, following the 
procedure of [^. We will not show the plots corresponding to model T-II here, since the 


13 















<—L 



5f 


>"T-< 

/ T™ \ 


j 


I ► ’ I ► 

(5) (5) 


# 




■ - X - 

Tx \ 


L 


L 


< I ► ■ I ► 

kQ Q k 

(H) {H) 


L 


FIG. 7: 2-loop neutrino mass diagrams in model T-II. The structure is as in BZ 


results are similar to the example models discussed in 
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From this exercise we estimate 


that typical uncertainties in the predicted proton-decay half-life should be around (2 — 2.5) 
orders of ma g^nit ude. The planned Hyper-Kamiokande experiment might cover a large part 
of this range 
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Note that, since model T-II has the same light leptoquark state T which 
appears in model T-I, the constraint on scalar mediated proton decay, see eq. ([9]), applies 
also in model T-II. 

We now briefly discuss neutrino masses. Since the GUT scale in model T-II is expected to 
be lower than in model T-I, the type-II seesaw contribution could be as large as ~ 6ylj^ 

meV, assuming /ta = ttla — rrtG- This is, in principle, large enough to explain the solar 
neutrino scale. However, 'mo and/or < 1 would render this contribution negligible. 

More important is the radiative contribution from TeV scale particles. The operators in 
eqs. (1T8|) and ffT9|) will produce 2-loop and 4-loop neutrino masses. The 

contribution to the neutrino mass matrix from 4-loop diagrams is expected to be at most 




0(10 


-lO'i 


eV 
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, so completely negligible numerically. The 2-loop diagram shown 
in £g. (I7j), on the other hand, will produce neutrino masses of the same order of eq. flTTll 
j^. We do not repeat the discussion of hts to neutrino data here, since it is very similar 
to that of the previous subsection. Also in model T-II LFV rates are expected to be near 
experimental limits. 

We turn now to the discussion of LHC searches. For the constraints on the leptoquark, 
T, see the discussion for model T-I. The same constraints apply, of course, also in model 


T-II. Cross sections for colour sextets are expected to be particularly large at the LHC [60 


since they can be produced as s-channel resonances. Both, ATLAS 6l|] and CMS 62| have 


searched for the appearance of new resonances in dijet spectra. Especially, 


62 


presents 


upper limits on a{pp —?■ jj) as function of the mass of a hypothetical resonance coupled to 
pairs of quarks. We can combine these limits with cross section calculations for resonances 
coupling to uu or dd pairs, see for example j3|. We estimate that, for my {mz) equal to 3 TeV 
this results in upper limits on the couplings (l/n^) of roughly < 0.19 < 0.056). 

Similar bounds can be derived for X (which nevertheless is an SU{2)l triplet, unlike Y and 
Z). Run-2 of the LHC will improve vastly on these numbers. 

Much smaller couplings (and larger masses) can be probed at the LHC, if there is a 
certain hierarchy in masses among the states X, Y, Z, R and T. We will discuss one 
example. If mx > 2mT, then the diquark X will have a non-negligible decay rate to a pair 
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of leptoquarks — see the top-right diagram in hg. ([S]). The branching ratio Br(X —)■ TT) 
depends on the relative size of the couplings h^^'^/rriT and y^j\ If the branching ratio is 
non-zero, since each of the T’s will decay into l~^j, the total signal is which enjoys 

a background orders of magnitude lower than the dijet spectrum. Since there is little or no 
background in this search for large values of mx, it might be possible to establish discovery 
with as few a (3-5) events. 

Given that we have no theory for the couplings, /rriT could be much smaller than 


f 7) 

Ulj . However, double beta decay depends on the same couplings. Thus, a measured hnite 
half-life in a future decay experiment would dehne a lower limit on as a function 

of the mass Mgff in a parameter range that should be completely coverable in the run-II of 
the LHC, considering the large diquark cross sections. 

Sensitivity for the doubly charged scalars at the LHC is much weaker, because of their 
colour singlet nature. Results of searches for same-sign dilepton pairs in electron and muon 


hnal states have been published by ATLAS 
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No deviation from SM expectations have 


been found and lower limits on the mass of R of roughly rriR > 400 GeV are derived for R 
decaying to either e^e^, e^/i^ or 
has a particularly rich phenomenology at the LHC. 


±^± ^±,,± ,,±,,± 100% branching ratio. In summary, model T-H 


IV. CONCLUSIONS 


li B — L violation exists at the TeV scale, one expects also sizeable contributions to the 
short-range part of the OudlS decay amplitude @, [^, with interesting consequences for the 


viability of leptogenesis |64l. I65l | . In this paper we have studied top-down scenarios inspired 
by SU{5) unihcation. These models violate B — L, thus generate neutrino masses and 
neutrino oscillation data can be easily htted. As in all TeV-scale models of loop neutrino 
masses, one expects that charged lepton flavour violation is large and possibly within reach 
of near future experiments. The presence of new states at the TeV scale changes the running 
of the three gauge couplings in such a way that, unlike in the standard model, gauge coupling 
unihcation can be achieved without introducing supersymmetry. 

We have constructed two exemplary models, one for each of the two tree-level topologies 
of 0i'(3(3 decay. Model T-I has a dark matter candidate and few low-scale (TeV) particles. 
On the other hand, there is model T-H which introduces several new scalars and thus has 
a richer phenomenology. The two models differ in their predicted gauge mediated proton 
decay half-lives, due to their different GUT scales. Model T-H might be in reach of the 
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planned Hyper-Kamiokande experiment 

The new states predicted by our models can be searched for at the LHC. We have discussed 
several different possible LHC searches, the most interesting of which are certainly the lepton 
number violating hnal states + 2j and + 4j, where I = e, fi, r. LHC can probe the 
part of parameter space of our models where the new states give the dominant contribution 
to the 0z//5/3 decay amplitude. In this sense, our models are falsihable experimentally. We 
are looking forward to the results of run-H at the LHC. 
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V. APPENDICES 

A — Lagrangians of the models 

This appendix contains the mass and interaction terms of the two models discussed in 
the main text. 


Model T-I 

Gonsider the usual three copies of the left-handed fermion representations Bp and lOi? 
of SU{5) in addition to one 24^. With an extended scalar sector consisting of one copy of 
the complex representations 5, 15, 45, 70 and a 24 real scalar field, the allowed mass and 


interaction terms are the following: 

■^mt = -^FF + -^FFS + -^SS + -^SSS + ‘^SSSS, (24) 

= ?T7,2424^24^-l-h.c., (25) 

^FFS = ^■%A0FJ5*+y^h0F^d0F,J5+^f5p,^10F,J45*+y^h0p,d0F,j^5 

+^f^5i.,i24j.70+^^°)24^24i.24+h.c., (26) 

^5 = m^5-5*+m^5l5-15*+m^545-45*+m2o70-70* + ^m2424-24, (27) 


if5g5=(^h(i)5.5.15* + h(2)5.i5*.70+h(3)5.24-45* + h(^)5-24-70* + h(®)l5-45*-45* 

15 ■ 45* • 70* + 15 • 70* ■ 70* + h(®)24- 45- 70* +h.c.) 5* ■ 24 

+ h(^°)l5-15*-24+h(^^’“) [24-45-45*]„^^ 2 +^^^^’“^ [70-70*-24]„^^ 2+^^^^^24-24-24, 
_ ’ ’ (28) 

^ssss = 5 ■ 5 ■ 5* ■ 70* + 5 • 5 • 45* • 70* + A® 5 ■ 5* ■ 45 ■ 70* + [5 ■ 45 • 45* ■ 70*] ^^^_2 

+^( 5 ;a) [ 5 * . 45 . 45 . 70 *]A( 6 ;a) [45 . 45 . 45 * . 70*]^^^ ... 5 + (other terms) , (29) 

where the Yukawa couplings and are symmetric under an exchange of the (ij) 
indices. The notation [•••]a=i 2 .. used above to indicate the existence of multiple gauge 
invariant contractions of a given product of fields. Quartic scalar couplings are numerous, 
therefore we do not write them all down here. The ones shown are important because, in 
their absence, it is possible to have a viable dark model candidate (which is the K scalar 
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field defined in the following; see the main text). We also note in passing that, as usual 
in SU{5), the scalar 5 and 45 are necessary to generate the Yukawa couplings 
which in turn will allow the SM fermions to have realistic couplings to the Higgs held (which 
is a combination of the SU{2) doublets found inside the 5 and 45). 

The SU{5) gauge symmetry breaks down into SU{3) x SU{2) x U{1) once the scalar 24 
acquires a VEV in the singlet component (under the SM group). Using the same notation 
as in the main text, the helds which remain light — besides the SM fermions — are 


hf = |^1,2,0 e 5 and 45, O = (8,1,0) e 24^., (30) 

T= (^3,2,0 e 15, K= (^1,4,0 €70. (31) 

The most general form of the interaction Lagrangian (including mass terms) which one can 
form with these helds is as follows: 

■^mt = ■^FF + -^FFS +-^38 +■^SSS + '^SSSS, (32) 

^FP = TnoO01 (33) 

^FFS = yffQiU^jH+y^fQd^^H*+yifLie<^^H*+yifL,d'jT+yl^^Q,OT* + h.c., (34) 

^ss = mjjHH*+m‘^TT*+m%KK*, (35) 

<^SS5 = 0, (36) 

^ssss= + + + +h.c.) 


Model T-II 

This model contains, as usual, three copies of the left-handed fermion representations 
5f and lOi? of SU{5). In addition, there are complex scalars 5, 45, 50, two copies of the 
complex scalar 15, and one copy of the real scalars 24 and 75. With this held content, we 
can write the following terms: 

■^int=-^FFS + .^SS + ^SSS + .^SSSS, (38) 

^FFS = ^f5FA0FJ5*+^fl0FA0F,J5 + ^f5F^^^^0F,,45*+^fl0F^,10FJ45 

+ ^pF,^5F,jl5k+^fl0FA0Fj50+h.c., (39) 

.^5 = m25-5* + (m25)..15i-15*+m2545-45*+m2o50-50* + ^m2424-24+^m2575-75, 

(40) 

5-5-15* + hlfl5i - 15j - 50+15* - 50* - 50* + h^^'> 5 - 24 - 45* + hf ^ 15* - 45* - 45* 
+/)(6;a) [24 - 45 - 45*] ^^^^2 24 - 45 - 50* + 5 - 45* - 75 + 5 - 50* - 75 
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■ 50* ■ 75+h.c.) 5 ■ 5* ■ 24+h[fh5i ■ 15* ■ 24+h^^^'>24-24-24 

+ /i(^^^50-50*-24+/J(^^’“) [45-45* ■75 ]„^i 2+^^^®^50-50*-75 + /J(^^^24-24-75 
+/^(i 8 ) 24 • 75 ■ 75+h^^^^ 75 • 75 • 75, (41) 

^ 5 S 5 = [15r 15,-• 50-24]^^, 2 + 15i ■ 15, • 50-75 +15i ■ 50* • 50* ■ 24 

+ Af[15i■ 50* ■ 50* ■ 75]^^, 2 + ' 15* ■ 24-24] 3 + (o^her terms). (42) 

Under an exchange of the {i,j) flavour indices symmetric, 

while and are anti-symmetric. As explained in appendix El the quartic couplings 
shown above are necessary in order to suppress neutron-antineutron oscillations in model 
T-II. 

With the helds 


H 


^1,2,-^ G 5 and 45, 
i^=(l,l,-2), X = 



(43) 

(44) 


which are assumed to remain light helds after of SU{5) symmetry breaking (see appendix 
El), one can build the following terms (note that there are two light T’s but only one Z): 


<^int — -^FFS + -^SS + -^SSS + -^SSSS, 

■H’ffs = yffQiU‘H + yfQd‘H' + y[fUe‘H’ + yftFd‘Tt 

+ y\fd‘d‘Z + slfe'e'il + +0,0,X + ylfu^u^Y + h.c., 

+ (mU - TiT* + m|ZZ* + mlRR* + m\XX* + m^YY *, 
^ggg = h\fTiTjX + h^^^YZZ + h^^^R*Y*Z + h^^'> X* X* Z + h.c., 

■^ssss = {X^^'^HHXZ + \^^^HHX*Y* + Ag),“^ [HT*T*Tl]^^^^ + xfTiT^X*Z 
+\^^^R*ZZZ + X^^'^RYYZ + [XXI *2 + X^^'^RX*X*Y + h.c.) 


+ X^^^HHH*H* + Xf^^HH*T,T* + X^^^^HH*ZZ* + X^^‘^'> H H* RR* 
+ [M*XX*],^,^2 + A(^^)m*XX* + Ag“) [T,T,T*Tn,=i,..._4 

+ xf '"^ [TiT*ZZ*]^^^^2 + >H/^RR*TiT* + xf'"'^ [TiT*XX*]^^^^ 

+ [TiT*XX*]^^^ 2 + + A(2^^i^i^*ZZ* 

+ A(22;a) [XX*ZZ*]^^,^2,3 + RRR*R* 

+ X^'^^^RRTXX* + X^^^'^RR*YY* + [XXX*X*]„^,_..._5 

+ [XX*XX*]^^, 2,3 + ■ 


(45) 

(46) 

(47) 

(48) 


( 49 ) 


There are the following symmetries in the parameters: 

• Vif j Vif 1 yij\ Vif j 8 ind A-^^ are symmetric under an exchange i j. 
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* ^ijk and have a mixed symmetry under permutations of the three indices. What 

this means is that for an S '3 permutation a, is equal to a linear combination of 

Ag’^) and Agf. The same is true for Agf.g. 

• Two of the Agg“^ parameters are symmetric both for an exchange of the {i,j) indices 
as well as (/c,/); the remaining two Agg“^ are anti-symmetric if we exchange the {i,j) 
or (/c,/) indices. 


B — Fine tunings 


Some parameters which multiply gauge invariant products of Gsm representations have 
relations among themselves in an S'f/(5)-symmetric theory. For example, in minimal SU{5) 
it is well known that yi = yj simply because decomposes as [d^Qj + LiCj) H* plus 

interactions involving the scalar colour triplet, which is heavy. There is no need to consider 
here the overall normalization factor in the contraction of the SU{5) representations but, 
on the other hand, the precise way of contracting the Gsm representations is important 
therefore let us briefly state how SU{2) and SU{3) indices are contracted. 

Take hrst SU{2). In the case of two doublets, we assume 2 ■ 2' = Cab^a^b where e is the 
Levi-Civita tensor (primes are used to distinguish fields transforming in the same way). If 
there is a triplet, one can picture its three components 3a (their electric charges depend on 
the hypercharge) as forming a matrix 


[3] 


be 




(50) 


in which case 2 • 2 ' ■ 3 = 2 . 2 ;[ 3 )^ and 2 * • 2 " ■ 3 = 

For colour indices, 3 • 3' ■ 3" = eabc^a^'b^c- Expressions with sextets are more easily 
expressed if we write the six components 6 a as a symmetric matrix: 


[ 6 ] 


be 


( 6i 

\ ^^65 


63 \ 


—6 


>6 


/ 


be 


(51) 


With this notation, 3" ■ 3'' ■ 6 = (3")J3'*)il6]„, and 6 • 6' • 6" = 

Finally, for both SU{2) and SU{3) indices, R - R* = Ra ■ {R*)a representation R. 

Having said this, we can expand some important SU (5) interactions and express them in 
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terms of the Gsm fields:® 

5F,ilOF,j5* oc {diQj + Lje^) H'* + (heavy), (52) 

oc [d^Qj — H”* + (heavy), (53) 

lOiT^jlOFjS oc {ulQj + Qiufj H' + (heavy), (54) 

10F,il0Fj45 oc {u^Qj — Qiuf^H" + (heavy), (55) 

5F,i5Fjl5 oc {d'^Lj + Lid'^j) T — + (heavy), (56) 

IOf.jIOfjSO (xulu^jY + QiQjX — + (heavy), (57) 

1545^50 oc TiTjX - Z^Z^Y + (heavy), (58) 

[15il5j50(24)]^ oc TiTjX + AZiZjY + (heavy), (59) 

1545^50(75) oc TiTjX + Z^ZjY + (heavy), (60) 

15i50*50* ocZ,(i^*r* + Y*R*) - ZiX*X* + (heavy), (61) 

15i50*50* (24) oc Z, {R*Y* + Y*R*)- ZiX*X* + (heavy), (62) 

15i50*50*(75) (x2Zi{R*Y* + Y*R*) + ZiX*X* + (heavy). (63) 


(The physical Higgs boson H is assumed to be a combination of the doublets H' and H" 
contained in the 5 and 45, respectively.) Interactions involving heavy Gsm representations 
are not shown. 

Equations (]5^ - (|55l) are important for obtaining realistic SM fermions masses; the 
Yukawa interaction in (l56|) participates in both our models, while (jS^ is only present 
in model T-II. This latter model contains scalar trilinear interactions 15 • 15 • 50 and 
15 ■ 50* ■ 50* as well, which lead to both Oi^/3/3 and neutron-antineutron oscillations. To 
suppress the latter process one can tune the trilinear coupling of 15 ■ 15 ■ 50 with the 
quartic coupling of 15 ■ 15 ■ 50 ■ (24) given that these two field contractions have different 
relative group-theoretical coefficients for the TTX and ZZY interactions. This is shown in 
equations fl58|) and fl5^ . On the other hand, 15■ 50* • 50* and 15 • 50* • 50* • (24) do share the 
same factors, so one must use instead the VEV of the 75 to suppress neutron-antineutron 
oscillations (see equations fl^ -fl63 |) ). 


In both models presented in this paper, one must also make sure that it is possible, 
by a suitable tuning of parameters, to make the various Gsm representations light (i.e., 
significantly smaller than the GUT scale). We shall illustrate this here only for T G 15 
although we have checked that the same is possible for all particles in both models. Consider 
then for simplicity a single copy of the 15, which contains T = (3,2,|), Z = (6,1, —|) and 


® These Georgi-JarIskog-like factors can be computed in a systematic way with the 
SubgroupEmbeddingCoefficients function of the Susyno program 6g. The user needs to inter¬ 
pret/adapt the output with care, since the program assumes a sign and normalization convention for the 
contraction of SU{5), SU{2) and colour indices which very likely differs from the one preferred/intended 
by the user. 
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A = Once the 24 acquires a VEV, the trilinear coupling 15 ■ 15* ■ 24 will contribute 

differently to the masses of these three Gsm representations: 

15 ■ 15* ■ ( 24 ) oc 4TT* - ZZ* - 6 AA*. (64) 

As far as quartic couplings are concerned, ignoring the scalar 75 representation 
(which may also acquire a VEV), one must consider three independent quartic cou¬ 
plings [15 ■ 15* ■ 24 ■ 24]^^j^ 2 3 (confer with equation (1421) b It turns out that 
[15 • 15* ■ (24) ■ (24)]^^^ 2,3 0^ Ca, iTT* -FCa,2 'Z'Z* -FCa,3 AA* with the vectors (Ca,l,Ca, 2 ,Ca, 3 )^ 
for a = 1,2,3 being linearly independent so, considering that linear combinations of the three 
gauge invariant contractions [15 ■ 15* ■ 24 ■ 24]^ are obviously gauge invariant as well, we 
may go ahead and define them such that 


[15-15* ■ (24) ■ (24)]^ oc TT*, 

(65) 

[15-15* - (24) - (24)]2 oc ZZ*, 

( 66 ) 

[15-15* - (24) - (24)]3 oc AA*. 

(67) 

Then, at the 517(5)-breaking scale. 

+ 4h (24) -h , 

( 68 ) 

m\ = — h.(24) -|- 

(69) 

"^Si. 3,1 = "^15 - 6h(24) + 

(70) 


where m\^ is the mass term of the 15 and h is the (properly normalized) trilinear coupling 
of 15 • 15* • 24 (corresponding to what we called in equations fl28l) and in fj4Tl) ). 
As such, it is clear that with an appropriate tuning of parameters, it is possible to make T 
and Z light while keeping A heavy. 


C — Decomposition table of 517(5) representations 


For reference, we provide here a table with the decomposition of SU (5) representations 
into those of the SM group (this data coincides with the one given in 67||). 


[1] F. F. Deppisch, M. Hirsch, and H. Pas, J.Phys. G39, 124007 (2012), arXiv: 1208.0727. 

[2] H. Pas, M. Hirsch, H. Klapdor-Kleingrothaus, and S. Kovalenko, Phys.Lett. B498, 35 (2001), 
arXiv:hep-ph/0008182. 

[3] J. Helo, M. Hirsch, S. Kovalenko, and H. Pas, Phys.Rev. D88, 011901 (2013), arXiv:1303.0899. 
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Rep. 


Decomposition 


50 

70 

70’ 

75 


35 

40 
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24 


5 

10 

15 



3,1 -f) + ( 6,2 + (3,2-1) + (8,3,0) + ( 8 , 1 , 0 ) + ( 1 , 1 , 0 ) + ( 6 , 2 ,|) + (3,2,|) + (3,1,|) 


TABLE I: SU{3) x SU{2) x U{1) decomposition of some SU{5) representations. As usual, the 



correctly normalized hypercharge yisf is obtained by multiplying the values of y shown here by 
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